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We study a noise-induced bifurcation in the vicinity of the threshold by using a perturbative 
expansion of the order parameter, called the Poincare-Lindstedt expansion. Each term of this series 
becomes divergent in the long time limit if the power spectrum of the noise does not vanish at 
zero frequency. These divergencies have a physical consequence : they modify the scaling of all the 
moments of the order parameter near the threshold and lead to a multifractal behaviour. We derive 
this anomalous scaling behaviour analytically by a resummation of the Poincare-Lindstedt series 
and show that the usual, deterministic, scalings are recovered when the noise has a low frequency 
cut-off. Our analysis reconciles apparently contradictory results found in the literature. 
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I. INTRODUCTION 



oo ■ 

o: 
o. 

vN ' A random noise can strongly affect the qualitative behaviour of a nonlinear dynamical system by shifting the 

' ~— j , bifurcation threshold, by modifying the characteristic exponents or by inducing unexpected transitions [l|, l3|. A 
I ^ ' straightforward approach to study the effect of noise on a bifurcation diagram would consist in analyzing the nonlinear 
^^ ' Langevin equation that governs the system. However, the interplay between noise and nonlinearity results in subtle 
—^ ■ effects that make nonlinear stochastic differential equations hard to handle [3,lj|- Therefore, most of the works either 
[ consider the linearized stochastic system in the neighbourhood of a stable manifold (the nonlinearity is thus eliminated) 
or analyze Fokker-Planck type evolution equations for the Probability Distribution Function (PDF) (the noise is thus 
integrated out by mapping a stochastic ordinary differential equation into a deterministic partial differential equation 
r) in the phase space of the system). Any type of noise can be studied by the linearization technique, however the 
behaviour of the linearized dynamics is in general different from that of the real nonlinear system. On the other 
hand, Fokker-Planck equations fully take into account the nonlinear dynamics but are valid only for white noise (the 
various 'Effective Fokker-Planck Equations' that have been proposed to treat colored noise have a restricted range of 
validity) . The Fokker-Planck approach can not therefore be applied to a noise with an arbitrary spectrum. 

The study of a nonlinear system subject to a deterministic forcing is a classical problem in the theory of dynam- 
ical systems (including, for example, the phenomenon of parametric resonance of an oscillator driven by a periodic 
modulation). Many mathematical methods have been developed to analyze this question, amount them the Poincare- 
Xj Lindstedt expansion which is a systematic perturbative calculation free of secular divergences SIgI- It is natural to try 

^^ , to adapt the Poincare-Lindstedt expansion to the case of nonlinear system driven by a noise [7|, |8[ : such an expansion 
T^ ' is versatile enough to deal with a nonlinear system subject to a noise having an arbitrary spectrum, thus allowing one 

I * to study precisely the impact of the spectral properties of the noise on a stochastic bifurcation. 

^-^ . The aim of the present work is to apply the Poincare-Lindstedt technique to the study of a stochastic Hopf 

QQ ' bifurcation and to determine the scaling exponents in the vicinity of the bifurcation threshold. We shall show that 
f^ , the characteristics of the noise have a strong influence on the scaling behaviour. In particular the existence of 
ILJ ' low frequencies in the noise power spectrum results in multifractality. The perturbative expansion also explains 
. ^H , the crossover between normal scaling and anomalous scaling and allows to resolve some controversial claims in the 
^ ' literature. 

J-j \ The plan of this work is as follows: In section Ull we define the model under study, carry out the Poincare-Lindstedt 

. . . 1 expansion and find the conditions under which this expansion becomes divergent. In section lllli we extract the leading 



term of this expansion to all orders and show that after resummation the normal scaling behaviour is replaced by 
multifractality. Concluding remarks are given in the last section. 

II. PERTURBATIVE ANALYSIS OF A NOISY HOPF BIFURCATION 

The random dynamical system we study here has been widely discussed as a paradigm for a noise-induced bifurcation 
[il, Q and is gouverned by the following stochastic equation : 

i:=(e + AC(t))x-x2P+i, (1) 



where p is a strictly positive integer and the noise ^(i) is a Gaussian stationary random process with zero mean value 
and with a correlation function given by 

(mm) = -Dili -t'\). (2) 

The power spectrum of the noise is the Fourier transform of the correlation function 

V{lu) ^ dtexp{iLut){^{t)^{0)) ^ dtexp{iujt)V{\t\). (3) 



When the noise ^(<) is white, the correlation function I? is a Dirac delta function and the power spectrum is a constant. 
We also recall that, because of the Wiener-Khinchin theorem Q, the function 'D{ijj) is non-negative. 
Applying elementary dimensional analysis to equation ([T]), we obtain the following scaling relations: 



X 



ti/i2p)^ e^<-i/2^ e^A^^f-i. (4) 

The dimension of the noise ^ is so chosen as to render the power spectrum "D^lu) a dimcnsionless function. 

A. The white noise case 

When ^(f) is a Gaussian white noise, the stationary solution of the Fokker-Planck equation corresponding to 
equation ([1]) is given by 

and where F represents the Euler Gamma-function. The bifurcation threshold is given by e = 0; for e < 0, the 
solution (O is not normalizable : the stationary PDF is the Dirac distribution 6{x) localized at the absorbing fixed 
point X = 0. For e > 0, the solution ([SJ is normalizable and is the extended stationary PDF. In this case, the moments 
of X are given by 

V{a) 
In the vicinity of the threshold, e is small and we find that the moments scale linearly with e, i.e., 

(a;2") ~ e {p^^y'/P-^V{n/p) . (7) 

In the following, we shall find the characteristics of the noise for which such an anomalous scaling is valid. 

B. The threshold shift 

The presence of noise can modify the bifurcation threshold which is given by e = in the deterministic case. For 
the stochastic differential equation ([T]) subject to an arbitrary noise ^(t), the critical value ec(A) that separates a 
localised PDF from an extended PDF, is determined by the vanishing of the Lyapunov exponent associated with the 
fixed point a; = 0. We linearize equation ([1]) around x = Q, 

— = (6 + AC(i))fa, (8) 

and find the Lyapunov exponent A to be 

A^ lim-(ln(<5.T)) = (^)=e. (9) 

t^oo dr ox 

This result is valid regardless of the nature of the noise ^ (we have only used the fact that the mean value of f 
vanishes). The Lyapunov exponent vanishes when e = 0. Thus, for the first order equation ([1]), there is no threshold 
shift due to the presence of noise and we always have ec(A) — 0. 



C. The Poincare-Lindstedt expansion 

In order to to simplify the general treatment, we shall study the case of a cubic nonlinearity, i.e., we take p = 1 in 
equation ([T]). The Poincare-Lindstedt method [a,H0jSl consists in writing two expansions : one for the solution x{t) 
of equation ([1]) and another for the deviation from threshold (e — ec(A)), both as power series of a formal parameter 
A, 

x{t) = \xi{t) + \^X2{t) + \^x:i{t) + X^x^it) + . . . , (10) 

e-ec(A) == Aei + A^ea + A^ea + A''e4 + . . . . (11) 

By substituting the expansions pUl and fTTI) in equation ([1]) we observe that x{t) and (e — ec(A)) are, respectively, 
odd and even in A (because equation JT]) is antisymmetric under x -^ —x). The Poincare-Lindstedt expansion thus 
reduces to (using the fact that ec(A) — 0) 

x{t) = X{xi{t) + X^X3{t) + \'^x^{t) + ...) , (12) 

e = A^ (e2 -h A2e4 + A^ee + . . .) . (13) 

Substituting these formal expansions in equation ([T]) and identifying the terms order by order in A, we obtain a 
hierarchy of first order differential equations for the functions X2i+i{t) 

Cx2i+i{t) =V {xi,xz,.. .,X2t-i,e2,...,e2t) , (14) 

where 'P is a polynomial function and C a linear differential operator : 

£=^-AC(t). (15) 

The hierarchy (fT4|) is solved with the initial conditions 

a;i(0) = 1 , and X2i+i(0) = for alH > 1 . (16) 

These conditions imply that a;(0) = A : the formal parameter A is simply equal to the value of x at time t — and, 
therefore, has the dimensions 

X^t^/^. (17) 

Besides, the numbers Ci appear as parameters in equation (|14p and are determined recursively thanks to the solvability 
condition that we now derive. Let us call yi{t) the solution of the adjoint equation C^yi — 0, which is given by 

yi{t)^exp(-Aj ^{u)duY (18) 

Multiplying equation ([T]) by the function yi and taking average values, we obtain 

(yiCxi) = {eyixi - yul) . (19) 

Integrating the left hand side of this equation by parts and taking into account the fact that yi is in the kernel of the 
adjoint operator C* , we derive the following relation 

e=j^. (20) 

By virtue of this solvability condition, the hierachy of equations (|14p is defined without ambiguity. The functions Xi 
and the parameters e^ are determined recursively in a unique manner using the initial conditions p6p . 
Eliminating A from equations (|12p and (I13p leads to the required expansion of x{t) in terms of e. 

D. Calculation of the first terms in the expansion 

We now calculate the first terms of the Poincare-Lindstedt expansion by applying the procedure described above. 
It will be useful to introduce the following auxiliary random variable Bt : 



Bt^ [ ^{u)du. (21) 

Jo 



Because ^ is taken to be a Gaussian random process, Bt is also Gaussian. The lowest order terms in the Poincare- 
Lindstedt expansion are then given by 

xiit) = exp(ABt) , (22) 

t 

2/ 



xsit) = xi{t)U2t~ / xi{{u)duj , (23) 

and the parameters £2, £4 are given by 

£2 = (a;?(i)) = (exp(2AS0), (24) 

£4 = 3(2:1X3) -£2 ( — ) , (25) 

\xi / 

where the expectation value (.) is taken over all the possible histories between times and t. 

E. Behaviour of the moments 

We now determine the behaviour of the moments of x, such as (a;^"), in the vicinity of the threshold, i.e., when 
£ — > 0. In the Poincare-Lindstedt expansion, this corresponds to the formal parameter A converging to 0. In this 
limit, we obtain 

(x2") = A^" ((xf) + 2nA2(xf -12:3) + ...)^ A2"(exp (2nABt)> = A^" exp {2n^A^B^)) , (26) 

e = A2(£2 + A£4 + ...)==; A2(exp(2ABt)) =A2exp(2A2(s2)) . (27) 

The last equality is derived by using the fact that Bt is a Gaussian random process. Eliminating A from equations (|26|) 
and (|27p. we obtain 



{X 



2n\ 



£"exp(2(n2-n)A2(s2)) . (28) 



This equation predicts a normal scaling behaviour identical to that of the deterministic case. In fact, this scaling is 
merely a logical consequence of the formal structure of the Poincare-Lindstedt expansion : indeed, we have (x^") r^ A^" 
and £ '^ A'^ and therefore (x^") ~ £". However, the proportionality factor between (x^") and £" in equation ([251 can 
be divergent when t — > 00. Such a divergence can change the scaling behaviour of the moments in the large time limit. 

F. Importance of lo^v frequencies in the noise spectrum 



In order to determine the multiplicative factor of £" in equation (|28p , we must calculate the variance of the random 
variable Bt : 



t rt i-t i-t 



{Bt)^ / / {^(unv))dudv ^ I / V{\u^v\))dudv^ / " 1 cosct gHdo;^ ^^9) 



The last integral is well defined at w = (the time t introduces an effective low frequency cut-off for oj ~ 1/t)- 
The behaviour of {Bf) for t — > 00 depends on the behaviour of I'(tt') at w ^ 0. The following two cases must be 
distinguished : 

(i) The spectrum of the noise vanishes at low frequencies, i.e., I'(O) = 0. Because I'(w) is an even function of w, 
we suppose that 2?(a;) ~ ui^ for oj ^ (we disgard non-analytic behaviour of the power spectrum at the origin. Such 
non-analyticity would correspond to long tails in the correlation function of the noise) . 

(ii) The power spectrum of the noise is finite at w = 0, i.e., X'(O) > 0. 

In case (i) , the long time limit of equation ([29]) is readily derived and we obtain (by using the Riemann-Lebesgue 
lemma) 



2, / V{ujy\u: 



Bt) -> / ' ' when t ^ 00 . (30) 



The variance of Bt has a finite limit at large times and, therefore, the prefactor of e" in equation ([28|) converges to 
a finite number when t — > 00, for all n > 1. Thus, the Poincare-Lindstedt expansion, used at the first order, leads 



to a well-defined asymptotic behaviour for (x^") and allows us to recover the 'normal' scaling behaviour which was 
predicted in [7] for the random frequency oscillator. Higher order terms in the expansion do not affect the scaling 
exponents and modify the prefactors only. The next order term was studied in [8[. 

In case (ii) , the integral on the right hand side of equation (|29p diverges when i — > cx3 and its leading behaviour is 



t r°° 1-cosu -- /-u 



— oo 



(5?> = z/ ^T^^(T)d.^-2?(0)i. (31) 



The variance of Bt grows linearly with time in the long time limit. Thus, for n > 2, the coefficient of e" in equation 
grows exponentially with time. Such a divergence of a prefactor is an indication of anomalous scaling. This anomalous 
scaling was not found in 8] where the authors analyzed only the first few terms of the Poincare-Lindstedt expansion 
and supposed, in order to avoid divergencies, that the noise spectrum has a low frequency cut-off uj^ (i.e., 'D{uj) = 
for uj < Lu^). We will show in the next section that even when low frequencies are present, the Poincare-Lindstedt 
expansion can be used to extract sound results and to predict anomalous scaling. 

III. RESUMMATION, ANOMALOUS SCALING AND INTERMITTENCY 

When the power spectrum at w = is finite {i.e., 2?(0) > 0), all the coefficients of the Poincare-Lindstedt se- 
ries (|12ll3p blow up when time goes to infinity. This divergence, that appears even for the lowest order, seemingly 
implies that the expansion breaks down when t ^ oo. We now show that it is still possible to use the Poincare- 
Lindstedt expansion and to extract from it by rcsummation the multifractal behaviour found in [l2|. 

A. Resummation of the Poincare-Lindstedt series 



Without loss of generality, we suppose in the sequel that I'(O) = 1; this amounts simply to redefining A as A/W I?(0). 
We first analyse the lowest order terms in the Poincare-Lindstedt expansion (|13I24I [25|) of e. If we retain only the 
most divergent contribution when i ^- oo, we obtain 

/„ A 9 N exp(8A2t) 3 exp(18A2i) 

62C.exp(2A2i), e,^ ^^ , .^ c. - -^^^^ . (32) 



We now investigate the general structure of the Poincare-Lindstedt series ([T3| , retaining for each order only the most 
divergent term. From dimensional analysis (equations (H]) and ^T7\ ). the dimensionless variable e/A^ can be written 
in terms of the dimensionless expansion parameter A/A as follows 

oo / a2 \ '^ 

^=E(-1)'"'«M A^ cxp(2fc2A2i); (33) 

k=l ^ ^ 

using equation ([5^ . we find ai = l,a2 = 1/2, 03 = 3/32. In fact, the general formula for the k-ih. coefficient Uk can 
be obtained recursively using equation ()14p 

"'= = 4^-i(fc-l)! • (^'^ 

We emphasize that the series ([SS| is divergent: its radius of convergence is strictly 0. However, it is possible to make 
a resummation of this divergent expansion by using the following representation for the exponential term 

2,_ f+°° du^ , 2 



exp(fc^) = / ^ exp(-u^ + 2fcu) . (35) 

J-00 V^ 

Inserting this identity in equation (133p . we obtain 

I^J -^exp (^-^) f:(-l)'=-^a, i^^y eM2kAVtu) . (36) 



A2 



k=l 



The series under the integral sign has generically a non-zero radius of convergence that depends on the coefficients 
afe. Defining 

oo 

F(z) = ^(-l)'=-iafcz\ (37) 

k=l 

and substituting in this equation the expression (|34p for the a^'s, we obtain 

F{4:z) = 4z(l - z) exp(-z) . (38) 

Equation (|36p can now be written as follows 

' ^^°° ^"^ :exp (-^] F f^cM'^AVtu)) . (39) 



A2 J_^ y/^K 

A resummation can be performed along the same lines for the mean value of a;^" : starting from the Poincare- 
Lindstedt expansion (fT^ . we find that 

/„2n\ °° , , / \\ 2("+fc) 

[X ; v-^ / , -.h, (n) -^ 



E(-1M"MX exp(2(n + fc)Wt). (40) 



fc=0 

Using equation ([55]) this divergent series is transformed as 



/2tt 
where we have introduced the new function 



exp -— G„ -^ exp(2AVt7.) , (41) 



Gn{z)^Y.^~l)H^^^z^+\ (42) 



__^^fe^,(")^«+fc 

fc=0 

To summarize, we have rcsummcd the Poincarc-Lindstedt series and obtained the following relations : 

— 9 ^ -^ ( A2' Ayi ) where T{z, w) ^ z when z — *■ and w is finite , (43) 



(x2" 



A2 

— Qn ( S2, A\/t J where t/„(2;, w) ~ z" when z ^ and w is finite . (44) 

The scaling behaviour of (x^") as a function of e in the vicinity of the bifurcation threshold is obtained by eliminating 
A between these two equations. 

B. Anomalous scaling 

The functions J- and Qn that appear in equations (|43l and I44|) have a singular behaviour when A ^ and t — > oo : 
these two limits do not commute. However, thanks to the expressions given in equations (|43[|44p . we can disentangle 
these two limits. 

If we take A ^ for a large but fixed value of i, we find from equations (j43p and (jH]) that (a;^") ~ e". Normal 
scaling is therefore recovered, in agreement with equation ([28]) . 

However, taking the limit t ^ oo first and keeping the value of A fixed and small, we obtain from equation 



e 



+ OC 



dv f v^ \^f^^ , •\ /■+°° Av ^/A2 , \ /+^dt;F(exp(i;)) 



exp -TTT^ F -^exp(^;) ^ / -— ==^ T7 exp(^;) - '^'°°,. ^ , (45) 



A2 J_^ 2Ax/2^ V 8A2t; V^^ / J-oo 2A^/2^ V^^ '7 2AV2^ 

where the last equality is obtained by the translation v ^r v — log(A2/A2). Using the explicit expression (|38p for F{z) 
we verify that the integral on the right-hand side of equation (|45p is a strictly positive real number, i.e., it is neither 
zero nor infinite. We have thus shown that 

' - '' (46) 



A2 AV2M 



where eo > . In a similar manner, using equation (|4ip. we find that the asymptotic behaviour of the 2n-th moment 
of X is given by 

/J2n\ 

(47) 




Ehminating t from equations (|46p and (|47p provides us the scahng behaviour of the moments of x in the vicinity of 
the bifurcation threshold 

(a;2") ~ C„eA2"-2 . (48) 

In the vicinity of e = all the moments scale linearly with e. This equation generalizes the calculation of the white 
noise case ([7]) to an arbitrary noise with non- vanishing zero-frequency power spectrum. We have thus shown, using the 
Poincare-Lindstedt expansion, that the low frequency of the noise strongly affect the scaling of the order parameter in 
the vicinity of the threshold of a stochastic bifurcation and induces a multifractal behaviour; a qualitative explanation 
of this effect was given in [I^l ■ 

C. Discussion of the random frequency oscillator 

The above analysis can be applied to more general random dynamical systems such as the parametrically driven 
damped anharmonic oscillator that naturally appears in the study of many instabilities [ij] . Such a system is described 
by the following equation : 

mx + m7i = (e + A^(t)) x ~ x^ , (49) 

where e is the control parameter and the modulation S,{t) is of arbitrary dynamics and statistics: it can be a periodic 
function or a random noise. For small driving amplitudes A, Liicke and Schank [7] have performed a Poincare- 
Lindstedt expansion, and obtained an expression for the threshold ec(A) (at first order in A). Their result has been 
verified both numerically and experimentally and is also in agreement with the exact result obtained for the Gaussian 
white noise (in this case a closed formula is available for ec(A) for arbitrary values of A). Another result obtained in 
P, la| is the scaling of the moments near the threshold, 

(x2") = sn (e - ec(A))" + O ((e - ee)"+^) , (50) 

where the constant s„ depends on ^[t) and on A. The moments have a normal scaling behaviour : (x^") scales as (a;^)". 
The bifurcation scaling exponent is equal to 1/2 and is the same as that of a deterministic Hopf bifurcation. However, 
this expression does not agree with the results for random iterated maps, for the random parametric oscillator and 
with recent studies on On-Off intermittency [10|, I 111 - Il2l | . These works predict that the variable x is intermittent and 
that the moments of x exhibit anomalous scaling, 

(x^") ~ K„(e - ec) for aU n > , (51) 

i.e., all the moments grow linearly with the distance from threshold. This multifractal behaviour, confirmed by 
numerical simulations for a Gaussian white noise, was derived using effective Fokker-Planck equations. 

The origin of the contradiction between equations (|50[) and (|5ip lies in the divergences that appear in the Poincare- 
Linsdtedt expansion. This fact, identified in [8], implies that the results of [7| are valid only for noises that do not have 
low frequencies. We have shown by studying a model technically simpler than equation (j49| . that the perturbative 
Poincare-Lindstedt expansion can be used for any kind of noise by resumming the divergent terms to all orders. This 
rcsummation allows to recover the scaling exponents for all the cases, and describes the crossover between the scalings 
given in equations (150p and ([5T|) . 

IV. CONCLUSION 

In this work, we have shown that the Poincare-Lindstedt expansion, a classical perturbative technique extensively 
used in the field of nonlinear dynamical systems, can be successfully adapted to analyze a stochastic model that plays 
the role of a paradigm for noise-induced bifurcations. These perturbative expansions for stochastic dynamics were 
studied in [7|, but, due to the presence of divergent terms, were applied only to random noises with vanishing power 
spectrum at low frequencies (e.g., noises with a low frequency cut-off) 8]. However, we have shown here that, by a 



resummation of the divergent terms, the Poincarc-Lindstedt method can be used for any type of noise. Moreover, the 
divergences that appear in the Poincare-Lindstedt expansion are not mathematical artefacts but have genuine physical 
consequences : the presence of low frequencies in the noise spectrum (that leads to these divergences) modifies the 
scaling behaviour of the order parameter in the vicinity of the bifurcation threshold. If low frequencies are absent 
from the noise spectrum, the order parameter has the same scaling as that of a deterministic bifurcation. In contrast, 
if the power spectrum does not vanish at zero frequency, the order parameter exhibits anomalous scaling in agreement 
with recent results for white and harmonic noise. 

Our work has allowed us to analyze precisely the role of low frequencies of the noise in a first order random dynamical 
system. The resummation technique we have used is fairly general and we believe that the results we have derived 
remain valid for systems of higher order in the vicinity of a stochastic Hopf bifurcation. 
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